arXiv:1509.01725v4 [math.DS] 1 Nov 2016 


On determinants of modified Bessel fnnctions and 
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Abstract 

We investigate the question on existence of entire solutions of well- 
known linear differential equations that are linearizations of nonlin¬ 
ear equations modeling the Josephson effect in superconductivity. We 
consider the modified Bessel functions Ij{x) of the first kind, which 
are Laurent series coefficients of the analytic function family e^O+l). 

For every I > 1 we study the family parametrized by k,n € Zb 
ki > ■■■ > ki, rii > ■■■ > Til oi {I X Z)-matrix functions formed 
by the modified Bessel functions of the first kind aij(x) = Ik^-mix), 
i,j = 1,... ,L We show that their determinants fk,n{x) are positive 
for every Z > 1, fc, n G Z* as above and x > 0. The above determi¬ 
nants are closely related to a sequence (indexed by 1) of families of 
double confluent Heun equations, which are linear second order differ¬ 
ential equations with two irregular singularities, at zero and at infinity. 
V.M.Buchstaber and S.I.Tertychnyi have constructed their holomor- 
phic solutions on C for an explicit class of parameter values and con¬ 
jectured that they do not exist for other parameter values. They have 
reduced their conjecture to the second conjecture saying that if an ap¬ 
propriate second similar equation has a polynomial solution, then the 
first one has no entire solution. They have proved the latter statement 
under the additional assumption (third conjecture) that fk,n{x) ^ 0 for 
fc = (Z,..., 1), n = (Z — 1,..., 0) and every x > 0. Our more general re¬ 
sult implies all the above conjectures, together with their corollary for 

‘Permanent address: Steklov Mathematical Institute, 8, Gubkina street, 119991, 
Moscow, Russia. Email: buchstab@mi.ras.ru 

^All-Russian Scientific Research Institute for Physical and Radio-Technical Measure¬ 
ments (VNIIFTRI), 

^Supported by part by RFBR grant 14-01-00506. 

SCNRS, France (UMR 5669 (UMPA, ENS de Lyon) and UMI 2615 (Lab. J.- 
V.Poncelet)), Lyon, France. Email: aglutsyu@ens-lyon.fr 

^National Research University Higher School of Economics (HSE), Moscow, Russia 
II Supported by part by RFBR grants 13-01-00969-a, 16-01-00748, 16-01-00766 and ANR 
grant ANR-13-JS01-0010. 


1 



the overdamped model of the Josephson junction in superconductivity: 
the description of adjacency points of phase-lock areas as solutions of 
explicit analytic equations. 


1 Introduction 

We consider the well-known problem on entire solutions of double conflu¬ 
ent Heun equations. Our results are directed to applications to nonlinear 
equations modeling the Josephson effect in superconductivity. 

1.1 Main result 

Let Y{l}) denote the space of the so-called two-sided Young diagrams of 
order 1: 

F(z') = {k = {h,... ,ki) \ h > ■ ■ ■ > ki} c . 

This notion is motivated by the fact that the k G Y{Zl‘) with fej > 0 are 
the classical Young diagrams. To every two-sided infinite number sequence 
and every / > 1 we associate the sequence of matrices nu¬ 
merated by pairs of two-sided Young diagrams k and n: 

Ak,n — ^ij — ^kj—riii fk,n ■— dfit Aj^n- (f'f) 

Remark 1.1 The matrices A^^n with k = (/ci, ki — 1,... ,ki — l + l) and n = 
(ni, ni — 1,..., ni — /-h 1) are the classical Hankel matrices [T5l pp. 301-302, 
495] written with inverse order of columns. The theory of Hankel matrices 
has important applications to the theory of functions, see 13 chapter XVI]. 
The matrices Ak^n corresponding to general two-sided Young diagrams may 
be considered as a natural generalization of Hankel matrices. The results of 
the paper provide the context on the crossing of dynamical systems, complex 
differential equations and physics where the latter matrices naturally arise. 

Remark 1.2 For every hxed two-sided Young diagram n the determinants 
fk,n with variable k G Y(Z^) form an infinite sequence of projective Pliicker 
coordinates corresponding to the subspace generated by the vector a = (aj) 
and its shifts by ni,..., in an infinite-dimensional space. 

The main result of the paper is the next theorem, which concerns the 
above determinants fk,n{x) constructed from the sequence of modified Bessel 
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functions of the first kind aj = Ij{x). Recall that Ij{x) are Laurent series 
coefficients for the family of analytic fnnctions 

+00 

gx{z) = ^ Ij{x)zK 

j=-oo 

Equivalently, they are defined by the integral formulas 

Ij{x) = — [ cos{j(j))d(j), j € Z. 

Jo 

Theorem 1.3 The determinant ^ix) in onp with aj = Ij{x) is positive 
for every two-sided Young diagrams k and n and every x > 0. 

Remark 1.4 Recall that a rectangular I x m-matrix is called strictly totally 
positive (nonnegative) [U fTHl ITfl [26] , if all its minors of all the dimensions 
are positive (nonnegative). Many results on characterization and properties 
of strictly totally positive matrices and their relations to other domains of 
mathematics, e.g., dynamical systems, mathematical physics, combinatorics, 
geometry and topology, are given in loc. cit. and in [271 [E] (see also 
references in all these papers and books). A.Postnikov’s paper |27j deals 
with the matrices I x m, m > I oi rank I satisfying a weaker condition of 
nonnegativity of just higher rank minors. One of its main resnlts provides an 
explicit combinatorial cell decomposition of the corresponding snbset in the 
Grassmanian G{l,m), called the totally nonnegative Grassmanian. The cells 
are coded by combinatorial types of appropriate planar networks. K.Talaska 
m obtained further development and generalization of Postnikov’s resnlt. 
In particular, for a given point of the totally nonnegative Grassmanian the 
results of m allow to decide what is its ambient cell and what are its 
affine coordinates in the cell. S.Fomin and A.Zelevinsky [12] studied a more 
general notion of total positivity (nonnegativity) for elements of a semisimple 
complex Lie group with a given double Bruhat cell decomposition. They 
have proved that the totally positive parts of the double Bruhat cells are 
bijectively parametrized by the product of the positive quadrant and 
the positive subgroup of the maximal torus. 

Theorem 11.31 provides an explicit one-dimensional family (given by clas¬ 
sical special functions) of strictly totally positive matrices Ak^n with I rows 
and infinite number of columns. We hope that appearance of strictly to¬ 
tally positive matrices in the context of the present paper would open a new 
direction of their applications. 
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Let us describe the relation of construction dLH) to Schur polynomials. 
Consider the generating function associated to a sequence aj: 

+ CXD 

M{a;w) = ajw\ rc € C. 

j = -OQ 

Recall that for every classical Young diagram n G Y {l}) one denotes 


One denotes 


^ j Zl') — 


^n2 


ni 




ni 


Z 

Z 


ni 

I 

n2 

I 


Z 


ni 

I 


I 

Mn{a-,z) = Aniz)Y\M{a]Zi), where z = (zi,... ,zi). 
i=l 


Lemma 1.5 One has the following formula 

Mn{a-z) = Y,fk,nZ^- ( 1 - 2 ) 

keT) 

Proof The Laurent coefficient at z^ of the function Mn{a\z) equals the 
sum (—taken over all the permutations a G Si. The 
latter sum is obviously equal to the determinant fk^n, see (11.11) . □ 


Remark 1.6 The determinants fk,n are well-defined for every k,n G . 
One has fk,n = 0 if either ki = kj, or = Uj for some i ^ j- If a 
tuple A: G is obtained from another tuple /c by a permutation a, then 
fkn ~ where sign((T) is the parity of the permutation a. 

Analogous statement holds for the other parameter n. The function Mn{a] z) 
is obviously anti-symmetric in z = (zi,..., z/). This together with Theorem 
11.31 implies that for every given two-sided Young diagram n and Oj = Ij{x) 
its Laurent coefficient with multi-index k is positive, whenever the order of 
the components kj differs from the decreasing one by an even permutation. 


Set 


5 = {I-1,1-2,. 


s\{z) 


^5{z) 


..,0) G Y(Z'), 

{zi,..., zi), A G Z^ 
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Recall that if A is a classical Young diagram, then by definition, s\ is the 
Schur polynomial associated to A, see [Ml P- 40]. If A G is not a Young 
diagram, then we will call s\ a Schur rational function. 

Corollary 1.7 For every A G Z^ one has 

Ms+x{a; z) = sx{z)Ms{a; z). 


Corollary 1.8 

^h-nia, z) — ^ ^ fk,n^kiz') — ^5(2;) ^ ^ fk^nSk—si^z) • 

k&Y{I}) fcgY(Zb 

The main application of our results concerns the family 

— sin (f> + B + A cos cut, w > 0, B > 0. 
at 


(1.3) 


of nonlinear equations, which arises in several models in physics, mechan¬ 
ics and geometry. For example, it describes the overdamped model of the 
Josephson junction (RSJ - model) in superconductivity (our main motiva¬ 
tion), see [Ml EQI EH m |28] ; it arises in planimeters, see [Ml Ell- Here u is 
a fixed constant, and {B, A) are the parameters. Set 

, B A 

T = ujt, 1 = F = IT- 

uj Zoj 

The variable change t ^ t transforms (|1.3I) to a non-autonomous ordinary 
differential equation on the two-torus x with coordinates i(p,T) G 

M2/27rZ2: 

• dd> sin N 

(j) = — = - \-1 + 2u, cos r. (1-4) 

or io 

The graphs of its solutions are the orbits of the vector field 


^ = — + I + 2lJ, cos T 

f = 1 


(1.5) 


on T^. The rotation number of its flow, see [21 p. 104], is a function p{B, A) 
of parameters. 

The phase-lock areas are the level subsets of the rotation number in the 
{B, H)-plane with non-empty interior. They have been studied by V.M.Buchstaber, 
O.V.Karpov, S.I.Tertychnyi et al, see mm and references therein and in 
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Section 3 below. Each phase-lock area is an infinite chain of adjacent do¬ 
mains separated by adjacency points. V.M.Buchstaber and S.I.Tertychnyi 
have described coordinates of a wide class of adjacency points [9] and con¬ 
jectured that this is the complete list of adjacencies. This was done via 
reduction of the family of non-linear equations to two families (13.41) . (13.9|) 
of second order linear differential equations of double confluent Heun type. 
They have shown that adjacencies correspond exactly to those parameter 
values, for which the linear equations have non-trivial holomorphic solu¬ 
tions at 0. They formulated conjectures on the existence of their holomor¬ 
phic solutions, which would imply the above conjecture on the adjacencies. 
It appears that the first linear equation cannot have polynomial solutions, 
while the second one can. The main conjecture, which implies the others, 
says that if the second linear equation has a polynomial solution, then the 
first linear equation does not have holomorphic solution at zero. Buch- 
staber and Tertychnyi have reduced it to their other conjecture saying that 
the determinants / 5 +(i,...,i), 5 (a:) are non-zero for every x > 0 [9l conjecture 
3, p.342]. They have proved their positivity for / < 3 and arbitrary x > 0 
and also for all I and small x > 0 dependently on I nm appendixes 1, 2]. 

Theorem Ol is a more general result, which implies all the above con¬ 
jectures. 

Theorem Ol will be proved in Section 2. Its proof is based on the 
following observation. The derivative of each modified Bessel function Ij is 
the half-sum of its neighbors Ij±i. This implies that the derivative of each 
determinant fk,n is a linear combination with positive coefficients of other 
determinants fk',n with k' obtained from k by adding ±1 to some component. 
The surprising idea is to write the latter formula for derivative not just for a 
single determinant fk,nj but for all the determinants fk^n with a fixed n and 
all k G Y{7}) together. This yields a linear ordinary differential equation on 
the infinite-dimensional vector function {fk,n{x))\k^Y('l}) with the right-hand 
side being equal to the discrete Laplacian plus the multiplication by 21. The 
latter right-hand side represents a linear bounded vector field on the space I 2 
of infinite sequences {fk)k€Y('i}) which the positive quadrant {fk > 0 | A: G 
F(Z0} is an invariant subset. It is shown that the infinite vector function 
{fk,nix))k£Y{i,‘) is / 2 -valued. The values at 0 of all the determinants being 
non-negative, the initial condition belongs to the positive quadrant. This 
implies that the above vector of determinants fk,n{x) lies in the positive 
quadrant for all x > 0, and hence, each determinant is non-negative for 
X > 0. Its strict positivity is deduced from the same differential equation, 
which implies that all its derivatives are non-negative and some of them does 
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not vanish at zero. 

Section 3 presents applications of Theorem 11.31 to double-confluent Heun 
equations and nonlinear equations (II.3p . (11.41) related to the Josephson ef¬ 
fect. It contains an introduction to the subject, an overview of previous 
results and proof of the above-mentioned conjectures of V.M.Buchstaber 
and S.I.Tertychnyi on Heun equations and adjacencies of phase-lock areas. 


2 Proof of Theorem 11.3 


In the proof of Theorem 11.31 we use the following classical properties of the 
modified Bessel functions Ij of the hrst kind, see [32l section 3.7]. 


Ij — I-j] ( 2 . 1 ) 

Ij\x>o > 0; lj{0) = 0 for j / 0; /o(0) > 0; (2.2) 

= /;. = i(/,_i + /,+i); (2.3) 

oo j-\-2s 

Ij (2y) = ^ - 0- (2-^) 


The next two propositions and corollary together imply that for every fixed 
n { 7 }) the vector function {fk,n{x))k£Y('Z‘) is a solution of a bounded lin¬ 
ear ordinary differential equation in the Hilbert space I 2 of infinite sequences 
{fk)keY{z‘)- ^ phase curve of a bounded linear vector field. We show that 
the positive quadrant {/^ > 0 | /c € Y{7^)} C I 2 is invariant under the 
positive flow of the latter field. This implies that fk,n{x) > 0 for all x > 0, 
and then we easily deduce that the latter inequality is strict for x > 0. This 
will prove Theorem 11.31 

Let us recall how the discrete Laplacian Adiscr acts on the space of 
functions / = f{k) in A: G For every j = 1,...,/ let Tj denote the 
corresponding shift operator: 


{Tjf){k) = fih, • • •, %-i, kj - 1, kj+i ,..., ki). 


Then 


Thus, one has 


i 


I^discr — + Hj ^ 

i=i 


2 ). 


(2.5) 


i 

i^^discrf^ip) ,Ps-l,Ps - ^,Ps+l,-- ■ ,Pl) 

S=1 
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+ f{pi, ■ ■ ■ ,Ps-l,Ps + l,Ps+l, ■ ■ ■ ,Pl)) - 2//(p). 


(2.6) 


Remark 2.1 We will deal with the class of sequences f{k) for which f{k) = 

0 whenever ki = kj for some i ^ j- It includes f{k) = fk,n{x) with fixed 
n € and x € M, see Remark 11.61 In this case the discrete Laplacian 
is well-defined by the above formulas (12.51) . (|2.6p on the restrictions of the 
latter sequences f{k) to the two-sided Young diagrams k G Y{7}). 

Proposition 2.2 For every I > 1 andn € Y{7}) the vector function {f{x,k) = 
fk,n{x))k^Y(i}) satisfies the following linear differential equation: 

Of 

^ = ^discrf + 21 f. (2.7) 

Equation (|2.7I1 follows immediately from definition, equation (12.31) and 
Remarks 11.6112.11 


Remark 2.3 For every k € Y{7}) the k-ih. component of the right-hand 
side in (12.7p is a linear combination with strictly positive coefficients of the 
components f{x,k') with k' G Y(Z^) obtained from k = (ki,...,ki) by 
adding ±1 to some k^. This follows from (12.6p . (12.7h . 


Proposition 2.4 For every constant R > 1 and every j > R? one has 


R^ 

Ij{x)\ < — for every 0 < x < R. 


( 2 . 8 ) 


Remark 2.5 The sequence ^ is bounded, and it decreases in j > R. 

Proof of Proposition 1^741 Fix an arbitrary j > Rf. Let us estimate the 
terms of the series (12.41) . For every s > 0 and y G [0, -^] one has 


yj+‘is ^j + 2s J^j I 

— - < ^ -< -^- 

s!(j-|-s)! 2.?+2s(j _|_ g)! j\ 2i+2« 


(R^Y ^ 1 R^ 
\TJ 


This together with (j2.4p implies (12.81) . 


□ 


Corollary 2.6 For every I > 1, n G Yflf) andx > 0 one has {fk,nY))k&(i}) ^ 
l 2 - Moreover, there exists a function C{R) = Cn{R) > 0 in i? > 1 such that 

^ \fk,n{xY < C{R) for every 0 < x < R. (2.9) 

kGY{l}) 
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Proof Fix an n G Y{7}) and an i? > 1 + |n|max, l^^lmax = niaxj \ni\. Set 

|^|n,max — |^|max l^lmaxj K = K{n^R) = \k G P(Z ) | |A;|ji^max ^ R }■ 

It suffices to prove uniform boundedness in x G [0, ii] of sum ()2.9I) taken 
through all k ^ K, since the complement Y(Z^) \ K is finite. Set 


M = max Ii(x). 
jGZ, 0<x<R 


The number M is finite, by (12.Sp and Remark 12.51 For every k a K one has 


^Ifcln.max 

\fk,ni^)\ < jrr, -for every 0 < x < R. (2.10) 

Indeed, if A: G /f, then some column of the matrix consists of functions 
Ij with j > |/i;|n,max > R^ , which satisfy inequality (j2.8p . by Proposition 
Ea and m- For the latter j the right-hand sides of inequality (12.8p are 
no greater than whenever x G [0,R], by Remark 12.51 The other 

matrix elements are no greater that M on [0, R]. Therefore, the module 
\fk,n{x)\ of its determinant defined as sum of l\ products of functions Ij 
satisfies inequality (|2.10l) . This implies that the sum in (|2.9I) through k ^ K 
is no greater than 


C{R) = 

keK 


^|^|n,max 

(|^|n,max)- 


< -|-C)0. 


The corollary is proved. 


□ 


Definition 2.7 Let H be the closure of an open convex subset in a Banach 
space. For every x G dH consider the union of all the rays issued from x 
that intersect H in at least two distinct points (including x). The closure of 
the latter union of rays is a convex cone, which will be here referred to, as 
the generating con^K{x). 

Proposition 2.8 Let H be a Banach space, LI C H be as above. Let v he a 
vector field on a neighborhood of the set Ll in H such that v{x) G K{x) 
for every x G dLl. Then the set H is invariant under the flow of the field v: 
each positive semitrajectory starting at LI is contained in LI. 

^The authors believe that this definition and the next proposition are well-known to 
specialists, but they have not found them in literature. 
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Proof Fix a point O G Int{Q). Consider the “anti-Euler” vector field w 
on H: its vector at a point x € H is the vector xO directed to O. Consider 
the family of vector fields = v + ew. For every e > 0 and x G dQ 
one has u^lx) G Int{K{x)). In other words, the field with e > 0 is 
directed strictly inside the domain fl, and its trajectories obviously stay 
in n. Hence, the latter statement holds for e = 0, by the theorem on 
the existence and uniqueness of solutions of bounded ordinary differential 
equation and continuity in parameter m- This proves the proposition. □ 

Now let us prove Theorem 11.31 Fix an n G Y{7}). The right-hand side 
of differential equation (j2.7p is a bounded linear vector field on the Hilbert 
space I 2 of sequences {fk)k&Y{i})- We will denote the latter vector field by v. 
Let Q C I 2 denote the “positive quadrant” defined by the inequalities fk > 0. 
For every point x G di} the vector v{x) lies in its generating cone K(x): the 
components of the field v are non-negative on fl, by Hemark l2.3l The vector 
function {fk,n{x))k^Y(i}) in a: > 0 is an Z 2 -valued solution of the correspond¬ 
ing differential equation, by Corollary 12.61 One has {fk,n{^))k&Y(i}) ^ 

/fc,n(0) = 0 whenever A; 7 ^ n; /n,n(0) = 4(0) > 0, (2.11) 

which follows from (12.21) . This together with Proposition 12.81 implies that 

fk,n{x) > 0 for every k G Y{l}) and x > 0. (2-12) 

Now let us prove that the inequality is strict for all k G Y{l}) and x > 0. 
Indeed, let fp^n{xo) = 0 for some p = {pi,... ,pi) G Y{7}) and xq > 0. All 
the derivatives of the function fp^n are non-negative, by (|2.7p . Remark 12.31 
and (I2.12p . Therefore, fp^n = 0 on the segment [0, xq]. This together with 
(EZI), Remark 12.31 and (|2.12p implies that fp/^n = 0 on [0, xq] for every p' 
obtained from p by adding ±1 to some component. We then get by induction 
that /n,n(0) = 0, - a contradiction to (|2.1ip . The proof of Theorem 11.31 is 
complete. 

3 Applications to double confluent Heun equations 
and Josephson effect: entire solutions and adja¬ 
cencies 

Here we prove the conjectures of V.M.Buchstaber and S.I.Tertychnyi from 
[ 9 ] mentioned in the introduction. They concern the family of nonlinear 
equations (|1.3I) : 

^ ~ sincp + B + A cos cut, A,uj > 0, B > 0. (3-1) 
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We fix an arbitrary a; > 0 and consider family (IHTTD depending on two 
variable parameters {B,A). The variable change t = ut transforms (|3.ip to 
the differential equation (|1.4p on the two-torus = 5^ x 5^ with coordinates 
{4‘jt) £ K^/27rZ^. Its solutions are tangent to the vector field 


0 + 1 + 2^ cos r 

f = 1 


B A 

oj 2uj 


(3.2) 


on the torus. The rotation number of the equation (|3.ip is, by definition, the 
rotation number of the flow of the field (j3.2l) . see [21 p. 104]. It is a function 
p{B,A) of parameters. (Normalization convention: the rotation number of 
a usual circle rotation equals the rotation angle divided by 27r.) The i?-axis 
will be called the abscissa, and the ^-axis will be called the ordinate. 

Definition 3.1 (cf. |18l definition l.lj) The l-th phase-lock area is the level 
set {B,A) I p{B,A) = /} C provided it has a non-empty interior. 

Remark 3.2 : phase-lock areas and Arnold tongnes. The behavior 
of phase-lock areas for small A demonstrates the Arnold tongues effect [21 
p. 110]. The phase-lock areas are called “Arnold tongues” in [18[ definition 
1 . 1 ]. 

Recall that the rotation number of system (|3.ip has the physical meaning 
of the mean voltage over a long time interval. The segments in which the 
phase-lock areas intersect horizontal lines correspond to the Shapiro steps 
on the voltage-current characteristic. 

It has been shown earlier that 

- the phase-lock areas exist only for integer values of the rotation number 
(a “quantization effect” observed in |6] and later also proved in [23119]); 

- the boundary of each phase-lock area {p = 1} consists of two analytic 
curves, which are the graphs of two functions B = gi^±{A) (see [7]; this fact 
was later explained by A.V.Klimenko via symmetry, see 122]); 

- the latter functions have Bessel asymptotics (observed and proved on 
physical level in [29]) see also [23] chaptrer 5], [S] section 11.1], [5]; proved 
mathematically in [22]). 

- each phase-lock area is an infinite chain of bounded domains going to 
infinity in the vertical direction, each two subsequent domains are separated 
by one point, the separation points lying outside the horizontal R-axis are 
called the adjacency points (or briefly adjacencies)-, 


11 


- for every I € Z the /-th phase-lock area is symmetric to the —/-th one 
with respect to the A-axis (symmetry of equation (13.11) : the set of adjacencies 
of all the phase-lock areas is also symmetric). 

There is a conjecture saying that for every I G Z the adjacencies of 
the /-th phase-lock area lie on the vertical line B = luj, see Fig.l. This 
conjecture is an open problem that was supported numerically in |18lpl . It 
was rigorously shown in loc. cit. that for every adjacency {B,A) one has 
I = ^ € Z, I = p{B, A){mod2) and |/| < \p{B,A)\. 

Theorem 13.121 one of the main results of the section stated and proved in 
Subsection 3.2 describes the adjacencies lying on a given line B = lu, I € Z, 
/ > 0, as solutions of an explicit analytic equation. To prove the conjecture, 
one has to show that their rotation numbers are equal to 1. 

The proof of Theorem 13.121 is based on the results of V.M.Buchstaber 
and S.I.Tertychnyi [9], Theorem 11.31 and the next theorem relating family 
(j3.1j) to families of double confluent Heun type equations. To state it, let us 
recall the following constructions and notations from [9]. Set 



UJ 





(3.3) 


The adjacencies correspond to / G Z, and it suffices to describe only those 
with / > 0, by symmetry. Thus, without loss of generality, everywhere below 
we consider that / G Z, / > 0. 

To family of equations ()3.ip . V.M.Buchstaber and S.I.Tertychnyi have 
associated in [9] the family of second order differential equations 

z‘^E" + ((/ + l)z + /x(l - z‘^))E' + (A - p{l + l)z)E = 0, (3.4) 

which is equivalent to the family of double confluent Heun equations 

z'^v" + ((/ -|- l)z -I- p,{z‘^ + l))u' -I- (A -I- p^)v = 0, V = e~^^E, (3.5) 


cf [9l equations (2), (23)]. 


Theorem 3.3 For every cu > 0, / G Z, / > 0 a pair {B, A) with A 0, 
B = l(jj is an adjacency for family of equations if and only if the cor¬ 

responding equation with A, p, as in I13.3\) has a nontrivial holomorphic 
solution at 0. 

^The results of paper [ 18 ] concern a slightly different family of differential equations 
equivalent to dMl), namely, ^ = i/sin 2 ;-|-a-|-ssinT. It is obtained from dMj) by coordinate 
and parameter change iot, x = —0, n=^,a=^,s=^ 
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ffdj=D.7 


i: 

□ 

vl 

11 

li 

iq 

Ji 

■» 

-11 

C 

□ 

•*•1 

ri 

i; 

bi 



-1.5 -1 -0.5 □ C.S 1 1.5 


Figure 1: Phase-lock areas and their adjacencies for oj = 0.7. The abscissa 
is -B, the ordinate is A. Figure taken from [9l p. 331]. 


Remark 3.4 A solution of equation f|3.4p is holomorphic at 0, if and only 
if it is an entire function: holomorphic on C. An entire solution is uniquely 
defined up to multiplicative constant. See m lemma 3, statement 4], 

Theorem 13.31 was implicitly stated in [9l p. 332, paragraph 2], We give 
its proof in Subsection 3.2 for completeness of presentation. 

In Subsection 3.1 we describe completely those parameter values for 
which equation ()3.4p has a nontrivial entire solution: these parameter val¬ 
ues are solutions of equation (|3.8I) . We then deduce the description of the 
adjacencies in Subsection 3.2. 
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3.1 Entire solutions of double confluent Heun equations 

For every A, /r G M, /i > 0 and / G Z, / > 0 set 

°° / 1 A \ 

Rm = Mj, where Mj = I 3 ( 3 - 1 ) ) m = / + 1, / + 2,..., 

j=m \ / 

(3.6) 

where the multipliers with bigger indices are placed to the right from those 
with smaller indices. The well-definedness of the infinite products Rm was 
proved in [Dj lemma 15]. Set 

6(A,^) = (A R^)Ri+i(^ J). (3.7) 

The main result of the present subsection is the following theorem. 

Theorem 3.5 For every I > 0 equation has a nontrivial entire solu¬ 

tion, if and only if the corresponding parameters (A, p) satisfy the equation 

^z(A,//) = 0. (3.8) 

Theorem 13.51 answers positively conjecture 2 from [9l p. 332]. 

Remark 3.6 Equation (13.4|] in an entire function E{z) can be translated 
as a system of linear three-term recurrent relations on its Taylor coefficients, 
see [H p.338, formula (34)]. The above matrices Mj and function ^;(A,/x) 
introduced in p. 337 ] originate from studying the latter recurrence rela¬ 
tions. The function fi(X, p) is analytic, and its Taylor series converges on the 
whole complex plane C^. This follows from the convergence of the infinite 
matrix product Rm, see [HI lemma 15]. (Formally, the proof of convergence 
in loc. cit. is given for A, /r G M, but it remains valid for complex values 
without changes.) 

The proof of Theorem 13.51 is based on the idea from [9] to consider 
simultaneously the equation 


z‘^E'' + ((-I + l)z + p{l - z^))E' + {p{l - l)z + \)E = 0, (3.9) 

which is obtained from equation (13.41) by changing I to —1. 

Theorem 3.7 theorem 7] The statement of Theorem \S.,7\ holds for I = 0. 


14 




Theorem 3.8 |3 theorem 2] For every Z € Z, Z > 0, and every (A, fi) 
satisfying equation \3.^) the corresponding equation \3.4\) has a nontrivial 
entire solution. 


Theorem 3.9 |3 theorem 8, p. 353] Let for a given I € N and some X, p, 
the corresponding differential equation 113.9\) have no polynomial solutions. 
Then equation i3.^) is also a necessary condition for the existence of a non¬ 
trivial entire solution of equation O- 

Theorem 3.10 If for a given I > 1 equation ii3.9\) has a polynomial solu¬ 
tion, then the corresponding equation has no nontrivial entire solution. 

Theorem 13.101 solves positively conjecture 1 in [9l p. 332], 

Proof Theorem 2 from [TO] says that the statement of Theorem 13. 101 holds 
under the additional condition that the determinant A{x) = fk,n{x), k = 
(I,... ,1), n = (Z — 1,..., 0 ) from (11.111 with aj = Ij{x) is non-zero for every 
a: > 0. But the latter inequality follows immediately from Theorem 11.31 
This proves Theorem 13. 101 □ 

Remark 3.11 The proof of theorem 2 from m, which relates the existence 
of polynomial solution of equation (|3.9h to the above determinant A(x), is 
done as follows. If equation ()3.4I1 has an entire solution, then each its solution 
is a holomorphic function on C* = C \ {0} whose Laurent series contains 
no monomial z®, —I < s < — 1 , see m lemma 3, part ( 6 )]. Suppose, 
by contradiction, that the corresponding equation (|3.9I1 has a polynomial 
solution E. Then degE = Z — 1, by [H remark 3, p.973]. The transformation 

9{z) ^ 

sends solutions of equation (13.911 to solutions of equation (13.411 [TOl lemma 
3, part (5)]. Therefore, the function 

E{z) = (3.10) 

is a solution of equation (I3.4|l . The vector of its Laurent coefficients at 
powers z®, s = —Z,..., —1 should be equal to zero. On the other hand, it 
is obtained from the vector of coefficients of the polynomial E (written in 
appropriate order with appropriate signs) by multiplication by the matrix 
Ak^n from of modihed Bessel functions aj = Ij{2p) with the above 

k and n: A{2p) = detAk^n- This follows immediately from formula (I3.10p . 
Therefore, if A(2/r) 7 ^ 0, then the above Laurent coefficient vector is nonzero, 
- a contradiction. 
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Proof of Theorem 13.51 , Theorem l3.5l follows from Theorems 13.TI - fSTTOl □ 


3.2 Josephson effect: adjacencies of phase-lock areas 

Theorem 3.12 For every given w > 0 and / € Z, / > 0, set B = loj, a 
pair {B, A) G with A ^ 0 is an adjacency of the corresponding family of 
equations if and only if the corresponding parameters X, p, given by 

US. 3\) satisfy equation F3. <j^) 

Theorem 13.121 follows from Theorems 13.31 and 13.51 
Proof of Theorem 13.31 Set 

$ = e*'^, z = F^ = e^^\ 


The complexified equation ()3.ip is equivalent to the Riccati equation 

^ = z-\{lz + piz"^ + 1 ))$ - ^(^>2 - 1 )). 

dz ziu 

The latter is the projectivization of the following linear equation in vector 
function {u,v) with = -: 

i — 2iuiz'^ (3.Ill 

I u’ = z-‘^{-{lz +p{l +z^))u + 


This reduction to a system of linear equations was earlier obtained in slightly 
different terms in [H |T3l O |20] . It is easy to check that a function v{z) is the 
component of a solution of system (j3.1ip . if and only if it satisfies double 
confluent Heun equation (j3.5p . or equivalently, the function E{z) = e^^v{z) 
satisfies equation (j3.4p . System (j3.1ip has singularities only at zero and at 
infinity; both are irregular ones. 

Let us suppose that given I > 0 and (A, p) correspond to an adjacency. 
Then the corresponding linear system (j3.1ip (and hence, equation (13.4|] ) 
has trivial monodromy operator along a positive circuit around the origin. 
This follows from the proof of m lemma 3.3]: it was shown in loc. cit. 
that the monodromy matrix should be equal to diag(l,e^™) with a G Z, 
and hence, to the identity. This implies that each solution of equation (13.4h 
is holomorphic on C* = C \ {0}. Exactly one non-trivial solution (up to 
multiplicative constant) should be holomorphic at zero. This follows from 
the fact that the germ at 0 of system (13.lip should be analytically equivalent 
to its diagonal formal normal form |18l lemma 3.3] and from [181 proposition 
2.9]. The first part of Theorem 13.31 is proved. 
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Now let us prove the converse. Let equation (lOD have a nontrivial 
solution holomorphic at 0. Then all its solutions are holomorphic in C*, by 
[9l theorem 3]. This implies that equation (|3.4I1 (and hence, system (|3.11l) l 
has trivial monodromy. This together with |18l proposition 3.2] implies that 
the parameters under consideration correspond to an adjacency. Theorem 
13.31 is proved. The proof of Theorem 13.121 is complete. □ 
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